Diffuse PeV neutrino emission from Ultra-Luminous Infrared Galaxies 
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Ultra-luminous infrared galaxies (ULIRGs) are the most luminous and intense starburst galaxies 
in the Universe. Both their star-formation rate (SFR) and gas surface mass density are very high, 
implying a high supernovae rate and an efficient energy conversion of energetic protons. A small 
fraction of these supernovae is the so-called hypernovae with a typical kinetic energy ~ 10 52 erg and a 
shock velocity > 10 9 cm s _1 . The strong shocks driven by hypernovae are able to accelerate cosmic 
ray protons up to 10 17 eV. These energetic protons lose a good fraction of their energy through 
proton- proton collision when ejected into very dense interstellar medium, and as a result, produce 
high energy neutrinos (< 5 PeV). Recent deep infrared surveys provide solid constraints on the 
number density of ULIRGs across a wide redshift range z ^ 2.3, allowing us to derive the flux of 
diffuse neutrinos from hypernovae. We find that at PeV energies, the diffuse neutrinos contributed by 
ULIRGs are comparable with the atmosphere neutrinos with the flux of 2 x 10 _9 GeVcm _2 s _1 sr , 
by assuming the injected cosmic ray spectrum as dNL/de' p oc £p~ 2 . 

PACS numbers: 95.85. Ry, 95.85.Hp, 98.70.Sa 



I. INTRODUCTION 

The galactic supernova remnants (SNRs) are widely 
suggested to be the dominant source for the cosmic rays 
(CRs) at energies below the "knee" at ~ 3 x 10 15 eV, 
most probably through the diffusive shock acceleration 
mechanism [lj . Though the details are still to be figured 
out, it is generally believed that the maximum energy 
of CRs accelerated by SNRs depends on both the ve- 
locity and the kinetic energy of the supernova outflow. 
A small fraction of the supernovae has a typical kinetic 
energy Ek ~ 0.5 — 5 x 10 52 erg and a typical velocity 
V ~ 10 9 cm s-^Sk/lO 52 erg) 1 / 2 (M SN /10M Q )- 1 / 2 , both 
are substantially larger than that of normal supernova, 
where Msn is the rest mass of the SN ejecta. These pe- 
culiar supernovae, such as SN 1997ef, SN 1997dq, SN 
1998bw and SN 2002ap, have been called the hypernovae 
0-0] • The maximum energy of the protons accelerated 
at the shock front of a supernova expanding into the uni- 
form dense interstellar medium (ISM) can be estimated 
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where n is the number density of ISM [5]. Such a 
fact motivates some colleagues to suggest that hyper- 
nova remnants are the dominant source of cosmic rays 
above the knee and the cosmic ray spectrum/flux up to 
~ 10 18 — 10 19 eV may be accounted for as long as the 
variety of supernovae has been taken into account [f34l0j . 

Ultra-luminous infrared galaxies (ULIRGs), first dis- 
covered in large numbers by the Infrared Astronomical 
Satellite in 1983, are among the most luminous objects in 
the local universe with infrared luminosity Ls-iooo^m > 
lO 12 I/0 [lUdH- The large infrared luminosities are at- 
tributed to large amounts of dust, which absorb ultra- 
violet (UV) photons and re-radiate them in the infrared 
(IR). Comprehensive observations show that ULIRGs are 
powered mainly by a large population of hot young stars, 



i.e., a "starburst", with a significant fraction also con- 
taining an IR-luminous AGN [13|. Local ULIRGs are 
exclusively mergers of gas rich galaxies accompanied by 
concentrated dust-enshrouded starburst, with very high 
star formation rate (SFR) > 200M Q yr -1 . The super- 
nova rate is expected to be high and huge amounts of 
cosmic ray particles are accelerated. Moreover, the rela- 
tively small sizes ~ 1 kpc of ULIRGs mean that the ISM 
density is also quite high, i.e., < n >~ 10 3 — 10 4 cm~ 3 [l4j|. 
suggesting that the accelerated high energy protons pro- 
duced in ULIRGs have a great chance to interact with 
the interstellar medium nucleons, to produce pions, and 
decay into secondary electrons and positrons, 7-rays, and 
neutrinos. If the energy loss time of the protons through 
the proton-proton collisions is shorter than the starburst 
lifetime and the confinement time, the protons will lose 
most of their energy before escaping and produce inter- 
esting observational signals [ill [l6| ■ 

Previously the possible high energy neutrino emission 
from the starburst galaxies has been investigated by as- 
suming that the observed GeV photons are produced via 
the decay of pions [l?], Gil • The events are expected to 
be significantly enhanced if the nearby starburst galax- 
ies host some Gamma-ray Bursts (GRBs) [l9j]. In this 
work, we focus on the PeV neutrino emission produced 
via the interaction of the accelerated Cosmic Rays with 
the dense environment in the ULIRGs. 



II. THE HYPERNOVAE RATE IN THE ULIRGS 

Though ULIRGs are very rare in the local universe, 
they are vastly more numerous at high redshifts. The 
relative contribution of ULIRGs to the SFR density of the 
universe also increases with redshift, and may even be the 
dominant component at z ~ 2 [2(|. Thus it is essential 
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to study the role of ULIRGs across a wide redshift range 
in producing PeV neutrinos. Observational data shows 
that core collapse supernovae of ty pe SN-Ib/c contribute 
with 11% of the total SN rate [21(. The hypernovae rate 
to the normal Ib/c SN rate is 7% in the local universe 
. Therefore the ratio of the hypernovae rate to the 
supernovae rate can be estimated to be /hs — 0.01. 

Recent IR observations show that the SFR density for 
ULIRGs increases rapidly at z < 1, and stays as a con- 
stant at higher redshift 1.2 ^ z ^ 2.3, which can be 
approximated as [23[ 

, , _ J 1.9 X 10- 4 e 3 - 3z M Q yr- 1 Mpc- 3 , z < 1.2; 
Psfr(z) - | o.OlMoyr-iMpc" 3 , 1.2 < z < 2.3. 

n (1) 

The supernova rate is related to the SFR via [24[ 

R SN (z) = 1.2 x 10- 2 psFR.(z)/Af Q . (2) 

Therefore, the hypernovae rate is 

Rhn(z) = /hs-Rsn(z) (3) 
2.4 x 10- 8 e 3 ' 3z yr- 1 Mpc- 3 , z < 1.2; 



1.2 x 10- b yr Mpc~ 3 , 1.2 < z < 2.3. 

Taking into account cosmology, the hypernovae number 
occurring in ULIRGs with redshift ^ z per year is nu- 
merically calculated by 

N MZ )=^r z dz w*)^)' 

#o Jo V^A + (l + z) 3 r! m ' 



Gas surface density can be derived from their SFR den- 
sity based on an empirical relation that SFR surface den- 
sity scales as some positive power /3 ~ 1.4 of the local gas 
surface density, i.e., the Kennicutt- Schmidt law (27| . 



Ssfr = (2.5±0.7)xl0" 



lM pc- 



1.4±0.15 



M Q yr 1 kpc 
(6) 

ULIRGs have SFR > 2OOM yr- 1 [2j|. If we assume a 
half-light radius of ~ 1 kpc, it yields a SFR density 
> 32 Mq yr _1 kpc -2 , consistent with observations of 
ULIRGs at both low and high redshifts. According to 
the Kennicutt- Schmidt law, the derived gas surface den- 
sity is > 1.0 g cm~ 2 . Hence, the energy loss time for 
the known ULIRGs is much shorter than the starburst 
lifetime in ULIRGs, which ranges from 10 7 to 10 8 years 

Another important factor to determine the fraction 
of protons energy conversion is the magnetic confine- 
ment time of protons. Adopting the similar confine- 
ment time as the Galaxy, and simply considering that 
the confinement time depends on the proton's Larmor 
radius(oc e'/B), the confinement time of the proton with 
energy e' p can be estimated as [ItT ] 
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where D e {z) = -tjl dz- 



It is found in Thompson et al. (2006) [30] that the 
magnetic field strength of the starburst galaxies scales 



distance, c is the light speed. Hereafter, we adopt the 
Hubble constant as Hq = 70kms~ 1 Mpc~ 1 , the matter 
density as 51 m = 0.3, and the dark energy density as 
Qa = 0.7 in the flat universe. 



is the comoving linearly with E gas , i.e., B w 3^Gt 
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quently, the confinement time can be rewritten as 
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III. DIFFUSE PEV NEUTRINO EMISSION 
FROM ULIRGS. 

The protons accelerated by the hypernovae would lose 
energy into the 7-ray photons, electrons and positrons, 
and neutrinos, through proton-proton collisions when in- 
jected into the interstellar mediums. A part of the pro- 
tons energy will convert into neutrinos via the decay of 
charged pions, ir + — >• /i + + — s- e + + v e + P M + and 
-> ii~ + P M -> er + v e + + v^. 

The energy loss time of protons is ri oss = (O.Sno-ppc)" 1 , 
where the factor 0.5 is inelasticity [25j |. and a pp is the 
inelastic nuclear collision cross section, which is > 70mb 
for the protons at energies e' p > 10 PeV that is of our 
great interest. Introducing a parameter S gas = m p nl as 
the surface mass density of the gas, with I as the scale 
of the dense region in the galaxies, the energy loss time 
reads HI 



Tlo 



= 1.4 x 10 4 yr- 



J gas 



100pc \ 1.0 g cm - 



(5) 



The fraction of the energy that the protons lose into pions 
is /tt = 1 — exp(— T con f/rioss)) which is close to 1 as long 
as T con f ^ ti oss . As a result, the protons with energy 
e' p lose almost all of their energy via the proton-proton 
collision before escaping from the starburst galaxies as 
long as ti oss < T con f, which constrains the critical gas 
surface density S C rit as 



S gas > S crit = 0.17 g cm 



lOPeV 



1/3 



lOOpc 



2/3 



(9) 

Consequently, for the ULIRGs with the gas surface den- 
sity S gas > 1.0 g cm -2 , the CR protons with energy up 

to - 2.0 x 10 3 PeV ( iooW) will lose almost all their 

energy via interacting with the dense ISM. 

The charged pion, whose energy is = 0.2£p, will 
then decay to produce four leptons, which share the en- 
ergy equally. Therefore, the fraction of the protons en- 
ergy converted into neutrinos is e„ = 0.05 [3l[. The 
observed neutrinos in the energy range (s v ,i, £vfi) are 
produced by the protons in the energy range (£^1,^.2) 
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in the ULIRGs, where e' p l = (1 + z)e Vt t/e v and e' 2 = 
(1 + z)e v _ 2 /^u- The energy fraction of the protons pro- 
ducing the neutrinos with energy between e Vl i and e v> 2 
is (for a > 2.0) 



c dcc 



~-p,2 
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where we assume the spectrum of the ejected protons 

dN' 

as -j^r oc e'~ a , and £p !min ~ 2GeV is the mini- 
mum energy of the ejected protons in the rest frame, 
and we define a parameter independent on the red- 

shift, edec = ,2-a /2-e. ■ Adopting an efficiency factor 

£ P,max" £ p m i n 

rj = 0.05 — 0.15 for the conversion of ejecta kinetic energy 
into the relativistic CR protons energy (3, [r| , the total 
energy of the CR protons is Ecu — 77-Ehn. Hereafter, 
we take the typical kinetic energy of the hypernova as 
£"hn = 2 x 10 52 erg and 77 = 0.1. Adopting a ~ 2.1, 
for the neutrinos with energies 0.5 — 5 PeV, we have 
edoc — 0.07. The total energy of the observed neutri- 
nos from e^i to e Vj % produced by each hyperonva in the 
ULIRGs is estimated as 



E v « 4xl0 4S erg 



2 x 10 52 erg 0.1 0.07 



(Ju_\ a ~ l f l + Z 

V0.05/ V 3 



(11) 

The produced neutrinos have the similar spectrum as 
the ejected protons, i.e., the observed neutrinos spectrum 
N c e~ a 31], then the normalized coefficient of 



is 



dN„ 
de v 



the neutrino spectrum can be calculated via 



N c . = 



E u (2-a) 



= AE V , 
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where we define a parameter A = 2-a~ a 2- Q to simplify 
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the expression. Consequently, the diffuse PeV neutrino 
flux integrating from local to the high redshift z reads 



F v {e u ) 



ATTD c {zYR w (z)N c ei 
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c /"* , Rw(z) 
dz- 



H0J0 (1 + z) a y/n M (i + z) 3 + n A 

where the luminosity distance Z?l — (1 + z)D c , while for 
the specified case with a — 2, it reads 



F v {e v ) = A*e'^ ec T]e u EuN 



X—— / dz 



(14) 



Rhn(z) 



H J (i + z) 2 y/n M (i + z) 3 + n A 



and e>* = ln ^~ ln ^ 



with A* = -, —, (uiu . 1 — -, — -, , — - 

i>rLZ v ^—ln.e u ^\ dec ' nE p m ax~™p m iE 

In figure 1., we show the flux of the diffuse neutrinos, 
at the energy of IPeV, from ULIRGs, GRBs and AGNs, 

for the ejected protons spectrum ^f- oc e'^ 2 . The diffuse 



PeV neutrino flux contributed by ULIRGs from local to 
z = 2.3 is F v ,p e y ^ 2 x lO^GeVou^ 2 s^sr" 1 for a = 2, 
which is comparable with the flux of atmosphere PeV 
neutrinos, implying an un-ignorable contribution to the 
PeV neutrino flux. However, the detection prospect of 
PeV neutrinos from ULIRGs is not promising. 

We estimate the amount of the detection rate of neu- 
trinos in the energy range (e 

V. 1 ) I/. ,2) via 

N v = Ae dec r]e^ l E w \ A exp (s v )s~ a ds u (15) 



dz- 



Rw{z) 



Ho Jo (i + zY'^n M (i + z) 3 + n A 

where A Qxp (e v ), varying with the energy of neutrinos, 
is the exposure coefficient of the detector for the diffuse 
neutrinos, with the unit of cm 2 s sr for a > 2. While for 
a = 2 we have 



N v = A*e^ cc rje u EuN / A e ^-p{e v )e v de v (16) 



dz- 



Rhn(z) 



Ho Jo (1 + z)V^m(1 + ^) 3 + ^a 

Considering the most sensitive neutrino detector nowa- 
days, i.e., the completed 86 strings IceCube observatory, 
adopting the effective area varying with the energy of 
protons, equation (fT6"|) gives N v = 0.1 for one year obser- 
vation. 



IV. CONCLUSION AND DISCUSSION. 

ULIRGs are a group of galaxies with ultra-luminous 
Infrared emission (Lir, > 1O 12 L0) and a high SFR 
(> 200Af Q yr~ 1 ). Consequently, a high hypernovae rate, 
Run ~ lOOOyr^Gpc" 3 at redshift 1.2 < z < 2.3 is ex- 
pected. Since the hypernovae can drive energetic shocks 
and are able to accelerate protons to 10 17 eV, we propose 

that huge amounts of protons with spectrum ^f- oc e'^~ a 

are ejected into these ULIRGs. The observations indi- 
cate that ULIRGs have a very high gas surface density, 
therefore the protons is expected to lose most of their en- 
ergy through interacting with the dense ISM in ULIRGs 
before escaping, providing a un-negligible contribution to 
the PeV neutrinos flux (~2x 10" 9 GeV cm" 2 s^sr -1 ). 
Its flux comparing with that of the atmosphere neutrinos, 
the GRB neutrinos and the AGN neutrinos have been 
presented in Fig.l. The ULIRG neutrino component is 
likely characterized by a cutoff (or break) at ~ a few PeV 
since the hypernovae are likely only able to accelerate the 
CR protons up to ~ 100 PeV and the ULIRGs can not 
confine the protons with energy much larger than ^100 
PeV, either. Such a component may be detected in 20 
years by the IceCube full configuration. 

Since the diffuse neutrino emission from ULIRGs are 
expected to be much rarer (see Fig.l), finally we suggest 
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FIG. 1: The flux of the diffuse neutrino emission from 
ULIRGs (purple solid line), GRBs (green solid line [321]). 
AGNs (red dotted line [18], [33j), assuming that the spectrum 
of the ejected protons is dNp/ds' p oc e' p ~ 2 ■ The black, red, 
green and purple dash-dotted lines represent the Greisen- 
Zatsepin-Kuzmin (GZK) neutrinos [34l . |35| referring to the 
models in [36|] (among Faranoff-Riley type II galaxies, i.e., 
FRII), 37] (with the best parameters that fit the cosmic ray 
data), |38| ] and [39l |. respectively. The black thick solid line 
represents the sensitivity of IceCube 86 strings for 5 years. 
The atmospheric neutrinos are presented by the data with 
error bars, which is measured by IceCube [4(j. The two 
black dash-triple-dotted lines are the upper bound and lower 
bound of the atmosphere neutrinos extrapolating to the high 
energy. 



that the two PeV neutrino candidates reported by Ice- 
Cube Collaboration [4l[, if cosmological, may be from 
other energetic sources, such as AGNs, GRBs (see [42j . 
but see (32L 43]), and cosmogenic neutrinos (see but 
see [45|, |46j). The origins of the reported two PeV neu- 
trinos are highly controversial so far, we anticipate more 
observations from the IceCube to draw a firm conclusion 
in the future. 
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